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Abstract 

We discuss a new solution, admitting the existence of dS/± branes, in five-dimensional 
Brans-Dicke theory. It is shown that, due to a special form of a bulk scalar field potential, 
for certain values of the model parameters the effective cosmological constant can be made 
small on the brane, where the hierarchy problem of gravitational interaction is solved. We 
' also discuss new stabilization mechanism which is based on the use of auxiliary fields. 
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q ; 1 Introduction 

Brane world models and their phenomenology have been widely discussed in the last years. The 
most known models - the Arkani-Hamed, Dimopoulos and Dvali scenario PQ and the Randall- 
Sundrum model [2J, - provide elegant, although different, solutions to the hierarchy problem of 
gravitational interaction. It seems that the second model is more consistent, because it takes 
^! into account the proper gravitational field of the branes. Meanwhile it was shown that the four- 
dimensional effective theory on the branes in the Randall-Sundrum model contains a massless 
scalar field - the radion, which is a consequence of the fact that the distance between the branes 
is not fixed by the parameters of the model. The coupling constant of this field to matter on 
the negative tension brane, which is assumed to trap the Standard Model fields, appears to be 
very large, which contradicts experimental data even at the level of classical experiments [31 H]. 

This problem was solved by introducing an extra scalar field living in the bulk. The most 
consistent model was proposed in paper [5] , where exact solutions to equations of motion for the 
background metric and the scalar field were found. The size of the extra dimension is defined 
O . by the boundary conditions for the scalar field on the branes. 

The models discussed above assume the metric on the branes to be the flat Minkowski metric. 
• t-h \ At the same time it is evident that more realistic models should account for a cosmological 
evolution on the branes. This problem is widely discussed in scientific literature, see, for 
example, papers [HIE], reviews [3(8] and references therein. Quite an interesting class of the 
brane world models is the one describing background solutions with dS^ metric on the branes. 
There are examples of such solutions in a slightly modified Randall-Sundrum model (with non- 
equal brane tensions) [HI [TU] . as well as in models with additional matter on the branes and in 
the bulk [11] including scalar field living in the bulk [T2J [13] . The latter models are of particular 
interest, because an additional scalar field can fix the size of the extra dimension thus giving 
stabilized models. 

One of the standard scalar-tensor theories of gravity is the Brans-Dicke theory (see [T^IT5]). 
In the context of brane world models this theory was also discussed in the literature for the case 
of static [16] and time dependent [T71 HE] solutions, including background solutions with <iS 4 
metric on the branes. Although one can transform the theory to the Einstein frame, in which 
the scalar field minimally couples to gravity, the theory in the original form can provide elegant 
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and simple background solutions. Moreover, in principle it is possible that we live in the world, 
in which five-dimensional scalar field is non-minimally coupled to five-dimensional curvature 
(this is defined by the interaction of matter on the branes with gravity, i.e. by the metric which 
we are supposed to perceive). In this case it is more convenient to consider untransformed 
action in the Jordan frame, in which the scalar field non-minimally couples to gravity. 

As it was mentioned above, brane world models are of particular interest because they 
provide elegant solution to the hierarchy problem of gravitational interaction. Nevertheless, 
such models should also describe cosmological evolution at least on the late stages. In this paper 
we discuss a stabilized brane world model in five- dimensional Brans-Dicke theory admitting dS± 
branes. We also discuss the values of fundamental parameters, which can make the effective 
cosmological constant on the brane which is supposed to contain SM fields (and where the 
hierarchy problem of gravitational interaction is solved) very small. 



2 The model 

Let us consider gravity in a five-dimensional space-time E = M4 x S l /Z2, interacting with 
two branes and with the scalar field 0. Let us denote coordinates in E by {x M } = {t,x l ,y}, 
M = 0,1,2,3,4, where x° = t; {x 1 }, i = 1,2,3 are three-dimensional spatial coordinates and 
the coordinate y = x 4 , —L < y < L, corresponds to the extra dimension. The extra dimension 
forms the orbifold S 1 /Z%, which is a circle of diameter 2L/tt with the points y and — y identified. 
Correspondingly, the metric gMN and the scalar field <fi satisfy the orbifold symmetry conditions 

9^(x, -y) = g^x, y), g^(x, -y) = -g^x, y), (1) 
#44(2:, -y) = g44(x, y), <f>(x, -y) = 0(x, y), 

H — 0, 1, 2, 3. The branes are located at the fixed points of the orbifold y = y = L. 
The action of the model has the form 



S = I d^x I dy-J^g 



-^-g MN d M <pd N <p-V(<p) 



(2) 



g\i(4>)d x — I ^—g\2(4>)d x. 
y=0 Jy=L 



Here V(<j)) is the scalar field potential in five-dimensional space-time, Xi^i^) are scalar field 
potentials on the branes, u is the five-dimensional Brans-Dicke parameter (we suppose that 
co 3> 1), g^v denotes induced metric on the branes. The signature of the metric gMN is chosen 
to be (— , +, +, +, +). Subscripts 1 and 2 label the branes. We also note that the dimension of 
the field is [mass] 3 . 

We consider the following standard form of the background metric, which is often used in 
brane world models (see, for example, [5]) 

ds 2 = lMN dx M dx N = e~ 2A{y) (~dt 2 + a 2 {t)r] t] dx l dx j ) + dy 2 (3) 

with rjij = diag(l, 1, 1), and the following form of the background solution for the scalar field 

<f)(x,y) = (f)(y) (4) 

(the background solution is the solution corresponding to the vacuum states on the branes, i.e. 
when all the fields on the branes are in their vacuum states). 
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In paper [16] a model, admitting a simple background solution in the case a(t) = 1, was 
proposed. Indeed, with 

V{4>) = M> (5) 

the background solution for the metric and scalar field takes a simple form [16] : 

A(y) = k\y\-kL, 

<j> = v(e- A )^, (6) 

where L is the size of the extra dimension, which is defined by boundary conditions on the 
branes together with constant v, and k is given by 

k 2 = -A, (tJ + 1)2 r. (7) 
(3cu + 4) (4w + 5) v ; 

Parameter A (and consequently parameter k) characterizes the energy scale of five- dimensional 
gravity. 

It appears that it is possible to modify slightly the bulk potential (J5J to get dS± background 
metric on the branes. We consider ansatz flUJ), (j3J) for the vacuum background solution with 

a(t) = e m , (8) 

( j ) = v(e- A )^, (9) 

where H is the four-dimensional Hubble parameter on the branes, v is a constant. The Hub- 
ble parameter H and the constant v will acquire their values after solving the corresponding 
equations of motion. Note that the form of the background solution for the scalar field is the 
same as the one used in [TB] for the case a(t) = 1 (see equation (jBj)). Assumption (jUj) simplifies 
considerably equations of motion for the scalar field and the Einstein equations following from 
©• 

The Einstein equations for a general form of the metric and the equation of motion for the 
Brans-Dicke scalar field can be found, for example, in |T7] (these equations are obtained for 
= ip 2 1 (8tt) and Ai j2 (0) = 0, nevertheless, the case Ai j2 (0) ^ can be easily restored). These 
equations have quite a complicated form. But for our choice of the background solution <^fy, 
OH]) and © the corresponding equations take a simpler form: 



1. /!Z/-component 



n 2A , 6^ + 8^ . /2 12w 2 + 31a; + 20 



2. 44-component 



QH 2 e 2A + r_L_,4" _ A 1 1 = (10) 

w + 1 (w + 1) 2 

^(0)^^(0) A 2 « 
= — — + — 7— o(y) + ——S(y ~ L), 




12H 2 e2 A_ Al ^ + 31u + 20^ m 

(U + 1) 2 ^ ^ 
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3. equation for the field 



1 ^ e » + ^±8 vt ._ ytfl 12^ + 3^ + 20 = (12) 

where ' — 4-. 

dy 

We suppose that the bulk scalar field potential has the form 

W) = A0 + ^r, (13) 

where A < 0, /3 < 0. The physical motivation of such a choice of the potential is the following: 
the first term is responsible for the solution to the hierarchy problem of gravitational interaction 
on the brane at y — L, whereas the second term is responsible for the constant Hubble parameter 
on the branes (i.e. dS^ space-time on the branes). One can see that the second term in ({TBI) is 
fine-tuned because the five-dimensional Brans-Dicke parameter u is utilized in the definition of 
potential (|T3|) . This fine-tuning of the five- dimensional bulk scalar field potential is similar to 
the choice of cosmological constant in four-dimensional gravity to get dS± space-time. Note that 
the vacuum structure of the four- dimensional space-time on the branes is defined by the bulk 
scalar field potential, not by the potentials on the branes. We also suppose that parameters (3, 
A are in the TeV energy range and do not contain hierarchical difference. 

First, let us consider equations ffTUl) . ffTTl) and ffT2]) in the bulk. The corresponding solution 
to these equations of motion has the form 

A = - In (de kM + C 2 e~ k \ y \) (14) 

with k defined by equation (J7|). We also can get 

H 2 = _4A; 2 |^C 1 C 2 (15) 
3w + 3 v ; 



and 



which means that the four-dimensional Hubble parameter is expressed through the parameter 
of the bulk scalar field potential, the five- dimensional Brans-Dicke parameter u and the 
constant v, which will be defined below. 

Let us discuss the values of the constants Ci, C 2 . We are interested in the effective theory 
on the brane at y = L (it will be shown below that the hierarchy problem is solved on this 
brane). To this end we should take such C\, C 2 that make the four-dimensional coordinates on 
this brane Galilean (see [31 IH HH]), i.e. A\ y= i = 0. Thus, 

de kL + C 2 e' kL = 1. (17) 

Equations ( 1T5|) and f lT7|) define the values of Ci, C 2 (we suppose, that the Hubble parameter 
H <C k because k is in the TeV range, whereas H should correspond to the late time accelerated 
expansion) . 

4(3o; + A)k 2 
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c ' K - ifc+W' (19) 



which results in 



e -A K e (^-%l) + |^±lg sinh (fcL - k\y\) = (20) 

= e<'-« + (^±Msinh(*L-%|). 

Now let us turn to the boundary conditions on the branes. They can be easily obtained 
from f flUj) . ffT2|) by the standard procedure (see, for example, |5j) and have the form 

1 ~ + I A \y=+o ~ 1^=0(0) , W 

o 6 ^ + 8 „/| A 2(0) 



LO + 1 



-A 



\y=L-0 J—\<t>=(j>(L), 



- -^-1^(0), (22) 

6^+8 rfA 2 (0) , 

z r~r A \y=L-o — J7 — \4>=4>(L)- 

00 + 1 d(j) 

Note that there are no boundary conditions coming from equation (TTTT) because it does not 
contain terms ~ A" and ~ S(y). 

We will consider the following form of the potential on the branes 



A0 2 8 



A 1)2 (0) = ±4 v / 3^+4^/ - T -^ T - t^J + M*) (0 - 0i, 2 ) , (23) 



where F\^{x) are scalar fields and <\>\p are constants. The absence of the kinetic terms for 
the fields F\^[x) looks rather strange. Nevertheless one can recall that supersymmetry is 
based on the use of such " auxiliary" fields, which are necessary for reaching the closure of the 
supersymmetry algebra [20] • A simple example with the fields of such type in classical field 
theory can be also found in [2"U] . 

Of course, one can use a more conventional form of the stabilizing potentials, for example 



■ A0 2 8 4(oo + lV3w + 4/3 , 

V 4^+5 PQi 

+ gi 2 (0-0l) 2 , 



, A0 2 8 4(w + lV3w + 4/3 . . 



4^+5 PV ^2 



+ 92 O- 02) 2 



The terms ~ gf 2 usually are introduced to ensure the absence of tachyonic modes in the 
linearized theory [21]. But such form of the potentials appears to be highly fine-tuned. One 
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can see that there is strong relation between the terms even if we consider only one of the brane 
potentials ( ([MP or ( I25p ). The auxiliary fields add extra degrees of freedom, which make the 
inner fine-tuning of the terms in the brane potentials unnecessary. Nevertheless, there remains 
fine-tuning in (123]) - the first term in (123]) is defined by the form of the bulk scalar field potential 
and by the form of the background metric ([3]) (one can check it by straightforward calculations). 
Such fine-tuning is inherent to almost all brane world models with compact extra dimension 
and two branes (see, for example, [21 15]). 

Equations of motions for fields F 12 (x) give (which can be obtained by means of the standard 
variation procedure with respect to the fields F\p[x)) 

0|,=o = 0i, (26) 

4>\y=L = 02- (27) 

We see that these equations do not contain fields Fx >2 (x) itself. Equations (12"T]) appear to be 
satisfied automatically for the choice (]2"3"]) (one can check it using ( ITT]) . (TT5|) and f lTB"]) ). whereas 
equations (122]) define the background values of the fields Fi j2 (x) 



= 4(, + 1^3^+4/3 



4( W + 1)^3^+4/3 (2Q) 



4^+5 

(we see, that the fields F\^{x) appear in equation of motion for the scalar field 0). Thus, 
although fields F\^[x) do not contain kinetic terms and are not dynamical, they can be treated 
as normal fields and one can use the standard variation technique to obtain corresponding 
equations of motion [2"U] . 

Equations (j2"6"]) and ([27]) define the size of the extra dimension. Indeed, from (J5]), (T2TJ]) and 
( 127]) it follows that v — 02- Then, neglecting the contribution of the term proportional to 
H 2 /k 2 (because we suppose that H <C k), we get from ([9]) and ( 126]) (see [16] for details) 

...(»+ DfcfM (30) 



A- 



'2 



Note that the same relation was obtained for the simpler model discussed in [16J. The four- 
dimensional Planck mass on the brane at y = L can also be easily obtained. For our purposes the 
contribution proportional to H 2 /k 2 in <$i§, ( ]2"U]) can be also neglected, and in this approximation 
the corresponding Planck mass is (detailed derivation can be found in [T6] ) 

M 2 Pl = \Jje- 2A dy~^ k e 2kL . (31) 

Thus, if kL 35 and parameters 0i, k of the theory lie in the TeV range, the hierarchy 
problem on the brane at y = L is solved in the way analogous to that used in the original 
Randall-Sundrum model [2]. Thus, we get weak four-dimensional gravity on the brane at 
y = L with Mpi ~ 10 19 GeV, whereas five-dimensional gravity is characterized by the TeV 
energy scale. 
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The five- dimensional Hubble parameter on the brane is defined by equation (j!6p . which can 
be rewritten as 

H 2 = -P^iM- (32) 

If we suppose that /3 ~ — lTeV 6u)+8 , u ~ 45 and 02 ~ 10 Te^ 3 , then the Hubble parameter in 
f )32|) has the same small value as that in ordinary four- dimensional gravity defined by vacuum 
energy density 

p A ~ H 2 M 2 Pl ~ 10~ 47 GW 4 . (33) 

Such value of H can correspond to the late time accelerated expansion of the Universe if we 
suppose that our four- dimensional world lives on the brane at y = L. 

Note, that such a small value of the effective cosmological constant appears because of the 
large power ~ 2u in the denominator of the fine-tuned bulk potential (|T3l) . One can argue that 
such a way of obtaining a small effective cosmological constant is analogous to introducing its 
small value "by hand". It is not the case in the model discussed above. It is somewhat similar 
to the solution of the hierarchy problem of gravitational interaction in the Arkani-Hamed, 
Dimopoulos and Dvali scenario [lj, where four-dimensional Planck mass on the brane appears 
to be large due to a large number of extra dimensions. But it is necessary to note that the large 
value 10 2aJ+2 is not introduced to the model by hand. Indeed, if 02 ~ ^TeV 3 , then from fl32|) it 
follows that H 2 w 15TeV 2 , which is extremely large value in comparison with the present day 
value of the Hubble parameter. Thus, the effective four- dimensional Hubble parameter on the 
brane depends on the value of the scalar field on the brane. The value 02 = lOTeV 3 does 
not create a new hierarchy itself. 



3 Stability 

A consistent study of stability of our model, at least in the linear approximation, appears to be 
quite a complicated task, which goes beyond the scope of this paper. Indeed, one should derive 
linearized equations of motion above the background solution, isolate the physical degrees of 
freedom and find the mass spectra of the excitations. Nevertheless we can simplify the problem 
and consider the simpler case = 0, for which we can show that our method of fixing the size 
of the extra dimension by utilizing the auxiliary fields F\^{x) does not lead to instabilities at 
least in the simplest cases. For (3 = the potentials on the branes and the bulk potential take 
the form (compare with those used in [16]): 

V{4>) = A0, A 1)2 (0) = ±4v^^^|0 + F 1)2 (x) (0 - li2 ) . (34) 

Boundary conditions f l2Tj) and fl22|) lead to the following background values 

Ft = 0, (35) 
F 2 = 0, (36) 

which also follow from fl2gp, fl2HD for (3 = 0. 

Linearized gravity in five- dimensional Brans-Dicke stabilized brane world models was thor- 
oughly examined in |21j. In particular, the model with 

V{<t>) = A0, A 1>2 = ±4 v ^A A /g±|0 + %0 - v lt2 ) 2 (37) 
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was also considered. It was shown that this model is stable and does not contain the scalar 
zero mode. The only difference between (13~4"|) and (13"T|) is the stabilizing potentials on the 
branes (such potentials are often called stabilizing potentials because they define the values of 
the scalar field on the branes, which leads to fixation of the size of the extra dimension and 
thus to stabilized model). The results obtained in [21J show that from the four- dimensional 
point of view linearized gravity can be described by the tensor (massless and massive tensor 
gravitons) and scalar (massive scalar modes) physical degrees of freedom. The tensor sector 
does not depend on the form of the stabilizing potentials and, as it was shown in [2T], does not 
contain tachyons. But the scalar sector of the model changes under the change of the stabilizing 
potentials. We will use results of |21j to show that the new method of stabilization does not 
lead to any new unwanted consequences in the scalar sector. We will not present here the full 
set of linearized equations of motion because these equations are quite tedious. One can find 
detailed calculations in [2Tj . 

To start with, let us parameterize the metric and the scalar field as 

9mn{x,v) = 1mn{v) + h MN (x,y), (38) 
<f>(x,y) = Mv)+ffav)> ( 39 ) 
Fi,2(x) = Fl 2 + j 1>2 (x)=j 1>2 (x). (40) 

For consistency with the previous sections below we will write and <p(y) for the background 
solution instead of F° 2 and <fio(y). 

It was shown in [21] that the physical degrees of freedom of the scalar sector for any form 
of the bulk scalar field potential V(<p) can be completely described by the new field g: 

g(x,y) = e-*U 2/3 (m*,!/) + |^ 

It was also shown that the fluctuations / of the stabilizing scalar field <fi can be expressed in 
terms of the field g through the gauge condition [21] 



The corresponding equation of motion for the field g in the bulk looks like [21] 

/ 15/3 2i\' 2 p 2A (h 5 / 3 e 2A 

( 9 'SH " s (3 " + 4) w> 9 + -ih^ 3 = °- (42) 



The difference between the model examined in [2T] and the model under consideration is the 
brane scalar field potentials. Indeed, from equations (T2"oT) . (I2"7j) it follows that /| y= o = f\ y =L = 
and thus using (jUJ) we get boundary conditions 

g'\y=+o = 0, (43) 
g'\ y =L-o = 0. 

The latter conditions differ from those obtained in [21] for (13T1) . 

There are also "boundary conditions", following from the linearized equation of motion for 
the field / [21], which, for our choice of the brane scalar field potentials fl34j) . take the form 

hi?) = - 3J ^e 2A d^g\ y=+0 , (44) 
j 2 (x) = ^e 2A d»d»g\ y=L ^. 
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We see from (j44l) that the fields ji{x), ]2{x) are not dynamical and appear to be defined by the 
values of the field g on the branes. The system appears to be not overconstrained. 
To make the mode decomposition, we represent g(x,y) in a standard way |21j : 



CO 



g( x , V) = ( Pn(x)g n (y), rfd^dvip^x) = /J> 2 n (p n (x). 



71=1 



Now the spectrum of scalar modes is defined by equations 



9: 



/ ( f> 5 / 3 e 2A \' 2 
n 0' 2 J ~9 



(3w + 4) 




9n + 




^n9n 







(45) 



9n\y=+0 — 0, 
9n\y=L-0 = 0, 



(46) 



where g n (y) is the wave function of the four- dimensional mode with the mass /i n [21]. It is 
not difficult to show (multiplying (l4"2j) by g> n and integrating it in the limits < y < L) that 
yU 2 > 0. As for the zero mode with /i = 0, its wave function g = 0. Thus, there is no zero 
scalar mode, which is inherent to stabilized brane world models. 

We showed that the model with j3 = is stable at least under the small fluctuations of the 
fields and our method of fixing the size of the extra dimension does not lead to instabilities. 
Other properties of the spectrum, such as orthogonality of the wave functions g n , can be easily 
obtained using the results of [21] . As for the case ^ 0, for a small H corresponding to 
one also expects stability of the model. Indeed, if all parameters of the model with /3 = lie in 
the TeV range, the masses of the lowest modes in the four-dimensional effective theory on the 
brane are also expected to lie in the TeV range, as well as the mass gaps between the modes. 
The case /3 ^ 0, providing the present day value of H, in principle leads to the modification of 
the spectra of tensor and scalar modes, but this modification can be neglected because of the 
extremely small value of H in comparison with TeV energy scale. In other words, we expect 
practically the same tower of the modes (including the massless graviton) propagating in the 
dS± space-time instead of the flat Minkowsky space-time, and it does not pose any problems 
with stability (we could expect ghost modes if the masses of the tensor modes were smaller 
than 2H 2 , as it happens in four- dimensional massive gravity [22], but it is not the case under 
consideration). Locally we can even neglect the influence of the non-zero Hubble parameter 
(for example, when considering the collider phenomenology or Newtonian gravity). For these 
reasons the model proposed in this paper seems to be stable. 

4 Conclusion 

In this paper we discussed stabilized brane world model in five- dimensional Brans-Dicke theory. 
The choice of the bulk potential was motivated by the demand to have dS^ space-time on the 
branes in the vacuum state (i.e. when there are no any matter fields on the branes). We 
note, that our solution is stationary, because the background solution for the scalar field does 
not depend on time, the size of the extra dimension is fixed and the four- dimensional Hubble 
parameter is constant. Such situation is realized in the limit x° — > oo. Indeed, in this case the 
ordinary and dark matter average densities on the branes tend to zero, solutions for the metric 
and the scalar field tend to this background solution. 
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Nevertheless the appropriate form of the bulk and brane scalar field potential should be fine- 
tuned in order to get dS± space-time on the branes. As it was mentioned above, such fine-tuning 
is the price we have to pay in order to get solutions with desired properties. The fine-tuning of 
the bulk potential is in some sense analogous to the choice of cosmological constant in ordinary 
four- dimensional gravity in order to get maximally symmetric dS^ space-time, whereas the 
fine-tuning of the brane potentials is necessary for the self-consistency of the solution. 

The advantages of the model presented above are the following. 

1. The hierarchy of gravitational interaction is solved in the model in the way analogous to 
the one utilized in the original Randall- Sundrum model [2] (see equation (1311 ). 

2. The size of the extra dimension is fixed (see equation (130]) ). 

3. The small four-dimensional Hubble parameter is defined by the parameter <f) 2 and can 
account for the late time accelerated expansion on the brane. 

We note that the effective cosmological constant on the branes appears to be small because 
of the large power in the denominator of the second term in the bulk scalar field potential 
(1T3"]) . which is similar in some sense to the solution of the hierarchy problem of gravitational 
interaction in the Arkani-Hamed, Dimopoulos and Dvali scenario [Ij. As for the fine-tuning of 
the scalar field potentials, in a most stabilized five-dimensional brane world models one should 
use special forms of the potential to get stationary vacuum solution, and the five-dimensional 
Brans-Dicke theory is not an exception. 

We hope that the results presented in this paper can be interesting for a future investigations 
of brane world models. 
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